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The exciton spin dynamics and polarization properties of the related emission are investigated 
in colloidal CdSe/CdS dot-in-rod (DiR) and spherical core/shell nanocrystal (NC) ensembles by 
magneto-optical photoluminescence (PL) spectroscopy in magnetic helds up to 15 T. It is shown 
that the degree of circular polarization (DCP) of the exciton emission induced by the magnetic field 
is affected by the NC geometry as well as the exciton fine structure and can provide information on 
nanorod orientation. A theory to describe the circular and linear polarization properties of the NC 
emission in magnetic field is developed. It takes into account phonon mediated coupling between 
the exciton fine structure states as well as the dielectric enhancement effect resulting from the 
anisotropic shell of DiR NCs. This theoretical approach is used to model the experimental results 
and allows us to explain most of the measured features. The spin dynamics of the dark excitons 
is investigated in magnetic fields by time-resolved photoluminescence. The results highlight the 
importance of confined acoustic phonons in the spin relaxation of dark excitons. The bare core 
surface as well as the core/shell interface give rise to an efficient spin relaxation channel, while the 
surface of core/shell NCs seems to play only a minor role. 


I. INTRODUCTION 

The on-going research and interest in colloidal NCs 
over the past three decades is partly driven by the 
prospect of building low cost, high quality and highly ef¬ 
ficient lighting or light harvesting devices^— such as col¬ 
loidal NC based light emitting diodes^ or solar cells^^ 
and partly by the opportunities opened up due to the 
high expertise in synthesis methods, enabling one to pro¬ 
duce NCs of almost arbitrary size and shape. A vast 
improvement regarding the optical properties of NCs has 
been the introduction of a shell, resulting in effective pas¬ 
sivation of surface trap states of the NC core, thereby 
considerably improving the quantum yield— It has been 
also shown that the shell provides the possibility to ma¬ 
nipulate the confinement of electrons and holes sepa¬ 
rately enabling to tailor their wavefunction overlap and 
thereby to control their interaction! 1 1_1 — In that way fur¬ 
ther improvement of the optical properties of NCs can be 
obtained. 

At low temperatures the light emission from NCs 
is predominantly determined by the exciton fine struc¬ 
ture, which in wurtzite type NCs arises from the inter¬ 
play of the crystal lattice anisotropy, the NC shape and 
the electron-hole exchange interaction— Due to its im¬ 
portance for the radiative emission of NCs the exciton 
fine structure has been studied theoretically and exper¬ 
imentally for various NC materials, sizes and shapes, 
such as spherical core-only and core/shell NCs*3r3i 
nanorods^Sr— and dot-in-rods (DiRs)— 

In spherical NCs the most crucial feature of the exci¬ 
ton fine structure is that the ground exciton state has 
total angular momentum F = ±2 and is therefore opti¬ 
cally forbidden in the electric-dipole approximation, i.e. 


is a dark state. Optically active states with F = ±1 
are split from the dark state typically by a couple of 
meV to higher energies. This splitting is often referred 
to as the bright-dark energy splitting and its magnitude 
depends on the electron-hole exchange interaction — No¬ 
tably it has been shown, that by increasing the CdS shell 
thickness due to the electron leakage into the shell the 
electron-hole exchange interaction and thereby the mag¬ 
nitude of the bright-dark splitting can be controlled— 
In contrast, recent experiments on DiRs indicate, that 
the shape of the CdS shell does not significantly affect 
the level order in the exciton fine structure for the core 
sizes of interest here 31— While in bare anisotropic struc¬ 
tures, such as nanorods, a change in the order of the 
valence subbands may occur at a certain diameter, lead¬ 
ing to an exciton ground state with total angular mo¬ 
mentum of F = 0, 2 - i 21 i 25 the fine structure of DiRs ap¬ 
pears to be similar to the spherical case. The order 
of the fine structure states is crucial for the photolu¬ 
minescence (PL) polarization of the NC, especially at 
low temperatures, where the thermal energy is smaller 
than the fine structure splitting and mostly the lowest 
lying exciton states contribute to radiative recombina¬ 
tion. While for DiR NCs high degrees of linear polar¬ 
ization of the PL have been experimentally observed at 
elevated temperatures*^— corresponding studies at low 
temperatures are still missing. 

The recombination dynamics of NCs is controlled by 
transitions between exciton fine structure states involving 
energy and spin relaxation. In order to conserve energy 
and angular momentum, these transitions are accompa¬ 
nied by absorption or emission of phonons. However, 
since phonons alone do not flip spins additional mech¬ 
anisms for spin relaxation are required. Understanding 
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these mechanisms is important because it enables design 
of spin-preserving environments for future spin-based de¬ 
vices in information storage . 29 ^ 30 Colloidal NCs, due to 
their high versatility regarding synthesis and processing 
as well as their unique and widely tunable optical prop¬ 
erties are appealing candidates for such applications.—^ 
Thus, comprehensive knowledge about the linkage be¬ 
tween NC size, shape and material properties and exci- 
ton spin dynamics is highly desirable. In recent exper¬ 
iments, exciton spin relaxation and its dependence on 
the NC shape have been addressed and possible mech¬ 
anisms involved have been discussed—2r— However, in 
these studies the authors consider only exciton relax¬ 
ation for bright-to-bright and bright-to-dark states and 
the experiments have been conducted at elevated tem¬ 
peratures. Transitions between the dark exciton ground 
states, which become important at low temperatures and 
occur in the presence of external perturbations, such as 
an external magnetic field, have not been discussed so 
far. As pointed out in Ref. HI, the dark exciton due to 
its long lifetime is a promising candidate for long-lived 
spin storage applications. 

In this paper, we present a detailed study of exci¬ 
ton spin dynamics and PL polarization in CdSe/CdS 
core/shell NC ensembles subject to an external magnetic 
field. We report on how the spin dynamics between the 
dark exciton ground states split by the magnetic field 
is affected by the size and shape of the CdS shell. To 
this end we measure the time-resolved degree of circular 
polarization (DCP) as function of temperature and mag¬ 
netic field for NCs covering a wide range of shapes and 
sizes, including DiRs as well as spherical core-only and 
core/shell NCs. Interestingly we find that the core/shell 
interface strongly affects the spin relaxation rate, whereas 
the shell shape and size have only slight impact in this 
regard. These results highlight the crucial role of the 
core/shell interface for the optical properties of NCs, es¬ 
pecially for the relaxation and recombination dynamics of 
dark excitons. Moreover, our results suggest that the spin 
flip is accompanied by the emission of confined acoustic 
phonons and remarkably NCs of spherical shape show a 
behavior similar to DiRs in this regard. 

We also study the impact of the shell shape on the PL 
circular polarization induced by the magnetic field. On 
the one hand we show that by measuring the PL-DCP one 
can probe the orientation of DiR ensembles, which after 
drop-casting on a substrate can show macroscopic align¬ 
ment features in so-called ” coffee stain rings” j 39 i 40 On the 
other hand, we analyze the dependence of the PL-DCP 
on magnetic field for DiRs and spherical core/shell NCs. 
Our results indicate, that the saturation level of the DCP 
in DiRs depends on the aspect ratio of the DiR shell. We 
elaborate a theoretical approach for analyzing and pre¬ 
dicting the polarization properties of NCs in the presence 
of an external magnetic field, which includes a dielectric 
enhancement effect resulting from the anisotropic shell 
as well as phonon-mediated coupling of the dark exciton 
ground state to bright exciton states with radiative tran¬ 


sitions parallel and perpendicular to the hexagonal axis. 
Interestingly our model suggests that, although at high 
temperatures the PL of DiR may be polarized along the 
hexagonal axis, at low temperatures and in the presence 
of magnetic fields it can be polarized perpendicular to it. 

The paper is organized as follows: In section [TT] in¬ 
formation on the synthesis of the NCs is given and the 
experimental setup is described. In section um the ex¬ 
perimental results are presented and discussed. Here, we 
start with a general discussion about how the orientation 
of wurtzite type NCs affects the circular polarization in 
magnetic field and we present results of measurements, in 
which such orientation effects are observed. Subsequently 
we explain the experimental procedure for investigating 
the exciton spin dynamics and the equilibrium PL-DCP. 
Thereafter results on the spin dynamics and DCP are 
presented and discussed. In section [TV] the theory for de¬ 
scribing the PL circular and linear polarization of DiRs 
in external magnetic field is presented and compared to 
the experimental data. The paper concludes with a sum¬ 
mary and an outlook in section m 


II. EXPERIMENTALS 

The wurtzite CdSe cores and the CdSe/CdS DiR col¬ 
loidal NCs were prepared following the protocol described 
in Ref. HT]. In this method, the use of phosphonic acids 
for the seeded growth of the CdS shell results in an 
anisotropic growth of the shell around the spherical core. 
The spherical CdSe/CdS NCs were prepared starting 
from the same wurtzite CdSe cores by using the recently 
reported ’’flash” synthesis. 42 This method makes use of 
oleic acid instead of phosphonic acids, which promotes an 
isotropic growth of the CdS shell and results in a spheri¬ 
cal morphology of the core/shell NCs. Parameters of the 
studied samples are collected in the Tables G] and HD for 
the DiRs and spherical NCs, respectively. 

Figure [T| shows typical transmission electron micro¬ 
scope (TEM) images of CdSe/CdS DiRs and spherical 
core/shell NCs. The geometrical parameters are illus¬ 
trated by the top sketches. The sizes of the CdSe/CdS 
DiRs and spherical core/shell NCs were determined from 
such TEM images. The size of the bare CdSe cores was 
determined from the energy of the first peak in the ab¬ 
sorption spectrum, using published data for the size de¬ 
pendence of the optical transitions^ 

For magneto-optical photoluminescence measure¬ 
ments, a solution of NCs in toluene with a concentration 
of a few /jM was deposited by drop-casting onto a glass 
slide and afterwards dried. For low temperature exper¬ 
iments the samples were mounted into an optical cryo¬ 
stat with direct optical assess to the sample via windows. 
The temperature could be varied between T = 2.2 K and 
20 K. The cryostat is equipped with a superconducting 
solenoid for generating magnetic fields up to 15 T. The 
magnetic field, B : was applied in Faraday geometry, i.e., 
parallel to the light wave vector and perpendicular to the 
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TABLE I. Parameters of the studied CdSe/CdS DiR samples. E m is the energy of the PL maximum, P® at - saturation level 
of the DCP, estimated by fitting Eq. © to the experimental data. g ex - exciton g factor, R e - enhancement factor, and 
r = R e P 2 o/P 2 i - parameter defined in Eq. ©. Numbers in brackets are standard deviations. 


# 

Core diameter 
D (nm) 

Shell length 
L (nm) 

Shell width 
W (nm) 

Aspect ratio 
L/W 

Em (eV) 

psat 

c 

Qex 

Re 

r 

1 

2.5 

28.9(1.5) 

5.3 (0.5) 

5.5 

2.18 

-0.65 

2.7 

3.2 

0.03 

2 

2.5 

22(1) 

3.8(0.3) 

5.8 

2.23 

-0.55 

2.7 

3.3 

0.13 

3 

2.5 

22.8(1.8) 

3.4 (0.2) 

6.7 

2.26 

-0.49 

2.7 

3.4 

0.19 

4 

3.2 

60 (5) 

4.1 (0.3) 

14.6 

2.11 

-0.63 

1.4 

3.8 

0.05 

5 

3.2 

32(1) 

4.1 (0.3) 

7.3 

2.12 

-0.60 

2.0 

3.5 

0.08 

6 

3.7 

20.4(1.3) 

6(0.6) 

3.4 

2.03 

-0.57 

1.4 

2.8 

0.16 

7 

3.7 

26.5 (2.7) 

5(0.3) 

5.3 

2.04 

-0.63 

1.5 

3.2 

0.06 


TABLE II. Parameters of the studied spherical samples. 


# 

Core Diameter D (nm) 

Shell thickness S (nm) 

Em (eV) 

psat 

c 

Qex 

Re 

r 

R1 

3.1 

2 

2.06 

-0.74 

2.7 

1 

0.02 

R2 

2.7 

0 

2.32 

-0.67 

1.2 

1 

0.25 

R3 

2.7 

0.8 

2.18 

-0.74 

2.4 

1 

0.02 

R4 

2.7 

3.7 

2.07 

-0.67 

2.8 

1 

0.25 

R5 

2.7 

5.3 

2.01 

-0.60 

0.9 

1 

0.5 



FIG. 1. (Color online) (a) Transmission electron microscope 
image of sample #4. (b) Transmission electron microscope 
image of sample #R5. The sketches on top illustrate the 
geometry of a DiR and a spherical core/shell NC. 


sample substrate. The sample was mounted on a holder 
top attached to a nanopositioner system, allowing for pre¬ 
cise positioning in the three spatial directions. NC photo¬ 
luminescence was excited with a linearly-polarized pulsed 
diode laser (pulse duration of about 100 ps), operating 
at a wavelength of 405 nm (photon energy of 3.06 eV) at 
a repetition rate of 1 MHz. The photoluminescence was 
detected either with an avalanche photodiode (time reso¬ 
lution of 100 ps) or spectrally-resolved by a monochroma¬ 
tor connected to a nitrogen-cooled charge-coupled device 
camera. The overall time resolution of the setup was 
determined to be 250 ps. The circular polarization of 


the PL was analyzed using a quarter wave plate and a 
Glan-Thomson- prism. 


III. EXPERIMENTAL RESULTS 


A. Effect of NC orientation on DCP 


Figure ©b) shows polarization-resolved PL spectra of 
sample #1 at B = 0 T and 15 T. Due to the size distribu¬ 
tion in the NC ensemble the PL linewidth is broadened 
up to approximately 80 meV (i.e. 20 nm). While at zero 
magnetic field the PL is unpolarized, it exhibits a strong 
DCP of —0.5 (i.e. —50%) at 15 T. The degree of circular 
polarization P c is calculated as 


-Lt+ 

= I a+ +1 (J — ’ 


(1) 


where I<7+ and I (j _are the PL intensities of the right- and 

left- circular polarized components, respectively. The 
spectral dependence of the DCP is shown in Fig. ©a). 
It changes only slightly across the PL band, which may 
be due to variation of the NC parameters or contribu¬ 
tion of spectral diffusion. Here we will not go into the 
details of these variations so that we can use spectrally- 
integrated PL intensities for calculation of the DCP. The 
observed polarization is result of an exciton distribution 
among the Zeeman-split sublevels, as illustrated in the 
inset of Fig. ©b). In NCs with wurtzite crystal struc¬ 
ture, due to the internal crystal anisotropy the exciton 
spin is constrained to orientations along the hexagonal 
c-axis of the crystal^ As a result, the magnitude of the 
exciton Zeeman splitting depends on the angle 0 between 
the crystal c-axis and the direction of the magnetic field: 
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FIG. 2. (Color online) (b) Photoluminescence spectra of DiR. 
sample #1 at B = 0 T and 15 T. Red and blue curves are 
the a + and a~ polarized components, respectively. The black 
curve corresponds to the unpolarized transition, as illustrated 
in the inset. The upper panel (a) shows the spectral depen¬ 
dence of the DCP. The insert shows the optical transitions 
between an exciton state |X) that is split in magnetic field 
into two spin sublevels and the unexcited ground state of the 
crystal |G). 


A E = g ex [iB^ ■ c = g e x^BB COS0, where g ex is the ex¬ 
citon g factor and gs is the Bohr magneton. The inten¬ 
sities I a ± in Eq. m therefore also depend on the angle 
0. The probability for detecting circular polarized light 
at an angle 0 with respect to the magnetic field axis is 
proportional to (1 ± cos0) 2 4i Thus, for an ensemble of 
NCs the P C (B) is derived by integrating Eq. ([]]) over all 
possible angles 

p , B s = /oVt+Qc) - I a -(x)]f OI {x)dx 
Jq [I<j+ (x) + I„-{x)\f OT (x)dx 

Here, x = cos© and f 0 i(x) is a weighting function de¬ 
scribing the orientation distribution in the ensemble. 

Assuming a thermal distribution between the Zeeman- 
split sublevels in equilibrium, the calculated DCP is 

pe= fo 2xtanh(AE/2 k B T)f OI {x)dx 
fo( 1 + X 2 )f or (x)dx 

where ks is the Boltzmann constant. For an ensemble 
of randomly oriented NCs the angle 0 between the crys¬ 


tal c-axis and the magnetic field direction is equally dis¬ 
tributed over all angles between 0° and 90°, i.e. f ov {x) = 
1. In this case, the integration over all angles 0 leads 
to a theoretical limit for P® q (B) of 0.75 (i.e. 75 %). 
However, for NCs with an anisotropic shape, such as 
nanorods or DiRs the situation can be different. For DiR 
NCs it is known that the c-axis of the core is typically 
parallel to the long axis of the rodi2£ Several methods 
have turned out to be suitable to prepare nanorod en¬ 
sembles with macroscopic alignment A 1,45- — Magnetoop¬ 
tical spectroscopy can be used for reliable characteriza¬ 
tion of the nanorod orientation in the ensemble by mea¬ 
suring P® q (P), as one can see from Eq. ([3]). In the case 
when all nanorods are oriented perpendicular to the mag¬ 
netic field direction, i.e. 0 = 90°, the expected DCP at 
any magnetic field P® q (P) = 0. On the contrary, when 
all nanorods are oriented along the magnetic field, i.e. 
0 = 0°, the resulting P® q (P) is maximal and depends 
only on the ratio AE/2ksT. This case corresponds to 
for(x) = 5(x — 1). For fully spin polarized excitons 
P® q (P) = 1. 

As a proof of principle we used this method to probe 
the nanorod orientation in ensembles prepared by the 
drop-casting method. When a solution of highly concen¬ 
trated NCs is drop-casted on a substrate, its distribu¬ 
tion on the substrate is subject to the so called “coffee 
stain effect” 1 39 i 4 ° 1 49 Due to this effect the nanoparticles 
accumulate in the outermost part of the drying droplet, 
leading to the formation of a coffee stain ring, where the 
NCs are densely packed. As it was shown in Ref. [39j in 
these coffee stain rings nanorods can align macroscopi- 
cally. Namely, they tend to be oriented parallel to each 
other and parallel to the substrate, whereas in the regions 
outside the ring the nanorods show no specific order and 
their orientation is random. 

Figure [3j(a) shows a PL intensity scan of a coffee stain 
ring section for sample ^5 measured by a confocal laser 
scanning microscopy with a spatial resolution of 10 /im. 
The spectrally-integrated PL intensity is plotted using 
a logarithmic grey scale with the brightness increasing 
with intensity. The high intensity within the ring sup¬ 
ports the dense packing of nanorods in this area. In the 
inner part of the ring (grey area) the nanorod density is 
substantially smaller. Outside the ring the intensity is 
almost zero but still measurable (black area). To obtain 
more insight the PL intensity profile along the red line 
in Fig. [3fa) is shown by the solid circles in Fig. [31(b). 
The time-integrated DCP at B = 5 T for the scan across 
the ring is also shown in this figure by the open circles. 
It is remarkable, that the DCP value being —0.23 in¬ 
side the ring drops almost to zero on the ring. This de¬ 
crease clearly evidences the nanorod alignment parallel to 
the substrate surface resulting from the vanishing exci¬ 
ton Zeeman splitting for the magnetic field perpendicular 
to the c-axis of the majority of nanorods. Note that the 
time-integrated DCP in general is lower than P® 9 . How¬ 
ever, if the spin relaxation time is much shorter than 
the exciton recombination time this difference becomes 
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negligible ! 44 ! 50 It will be shown below that this is always 
the case for the samples investigated here. 

The technique presented is suitable to yield informa¬ 
tion about the nanorod orientation only with respect to 
the magnetic field axis perpendicular to the sample sub¬ 
strate. Although it cannot compete with electron mi¬ 
croscopy in terms of resolution, it has noteworthy ad¬ 
vantages. While for TEM imaging special carbon grids, 
which are transparent to the electron beams, have to be 
used, the method we propose is applicable to a wide range 
of substrate materials and can be used for device charac¬ 
terization. Furthermore, while the area available to TEM 
is typically limited to a few microns, the optical tech¬ 
nique can be applied to much larger scales. Its spatial 
resolution is limited to several hundred nanometers due 
to light diffraction, depending on the laser wavelength 
and the numerical aperture of the optical excitation and 
detection system. 


B. Experimental investigation of equilibrium DCP 
and exciton spin dynamics 

In the preceding section we have highlighted the im¬ 
pact of the nanorod orientation on the magnetic-field- 
induced DCP in DiR nanocrystals. Therefore, when in¬ 
vestigating the magneto-optical properties, this aspect 
has to be taken into account. In order to ensure com¬ 
parability between different samples, all measurements 
presented in the following were performed in the center 
of the dried droplets, where no specific order of the DiR 
particles is expected. 

The internal crystal structure and NC shape 
anisotropy, in combination with the electron-hole ex¬ 
change interaction, lift the eight-fold degeneracy of the 
exciton ground state, leading to an exciton fine structure 
consisting of five energy levels, which are distinguished 
by their total angular momentum^ At low temperatures 
it is typically sufficient to consider only the lowest exci¬ 
ton states when discussing the optical properties of NCs. 
It was shown in previous studies^ that the exciton fine 
structure of DiRs for the core sizes of interest here is quite 
similar to the case of spherical core/shell NCs. A sketch 
of the lowest fine structure states is shown in Fig. Oja) 
without and in presence of an external magnetic field. 
The lowest exciton ground state |F) has total angular 
momentum of |F| = 2 and is optically forbidden. It lies 
below an optically allowed state | A) with |F| = 1. The 
bright-dark energy splitting A Etd is strongly dependent 
on the shell thickness, because a change in shell thick¬ 
ness leads to a modified electron-hole exchange interac¬ 
tion. This aspect was investigated in several experimen¬ 
tal studies^ ! — i — For the samples we investigate here the 
magnitude of A Em lies in the range of 2—5 meVi 24 At low 
temperatures PL emission occurs mainly from the dark 
state, due to efficient spin flip scattering from the bright 
to the dark exciton state, which occurs on picosecond 
timescalei 18 i 51 Dark exciton radiative recombination is a 


long standing question and its origin has not been fully 
disclosed. Several reports evidenced a phonon-assisted 
mechanism involving longitudinal optical and longitu¬ 
dinal acoustic phonons. This, however, cannot explain 
some peculiar experimental observations, such as pres¬ 
ence of the zero phonon lin o 15 i 23 i 52 i 53 or shortening of 
the PL decay with increasing external magnetic fields 

In the presence of a magnetic field the exciton states 
are split into Zeeman doublets with well-defined spin ori¬ 
entation. At low temperatures the establishment of ther¬ 
mal equilibrium population among the exciton spin sub- 
levels involves the bright-dark spin relaxation character¬ 
ized by a time taf and the spin relaxation between the 
Zeeman sublevels of the bright and dark excitons. Since 
exciton relaxation between |+2) and |—2) requires a si¬ 
multaneous spin flip of electron and hole, it is expected 
that this process is less efficient compared to the bright 
to dark exciton relaxation, i.e. r s taf , where r s is 
the spin relaxation time of the dark exciton^ 7 The spin 
relaxation between the |+2) and |—2} states is the main 
focus of this study. 

The spin dynamics can be investigated by analyzing 
the time-resolved DCP P c (t ), as shown in Fig. SJb). Af¬ 
ter pulsed excitation, the a + and cr _ polarized compo¬ 
nents are traced as function of time. The time-resolved 
DCP, which is shown by the black curve, is calculated 
from the experimental data after Eq. ©■ Right after ex¬ 
citation the PL is unpolarized, because the photoexcited 
excitons being excited with linearly polarized (unpolar¬ 
ized) laser light with large excess energy do not have 
a preferential spin orientation. Subsequently the exci¬ 
tons relax into the energetically lower lying spin states 
until they achieve thermal equilibrium, as evidenced by 
the plateau of the DCP in Fig. ©b). The characteris¬ 
tic time of the spin relaxation process involved can be 
quantified by determining the characteristic rise time of 
the time-resolved DCP^i In the following we analyze the 
dependence of the PL polarization and spin dynamics on 
magnetic field, temperature and shape of the NCs. To 
that end, we fit the time-resolved DCP curves with the 
equation 

p c{t ) = ^ c eq [l - exp(— t/r s )\, (4) 

where P® q denotes the DCP in thermal equilibrium and 
t s the spin relaxation time of the dark exciton. A typical 
fit is shown by the green curve in Fig. ©b). 

Figures ©c) and©d) show the P c {t) at various B and 
T for sample #1. On the one hand, as the tempera¬ 
ture is kept constant and the magnetic field is increased, 
the DCP saturates at higher values, due to the increased 
magnitude of the Zeeman splitting. On the other hand, 
the spin relaxation process becomes more efficient, as can 
be seen from the faster DCP increase at higher fields. If 
the magnetic field strength is fixed and the temperature 
is increased, the saturation level decreases due to thermal 
population of the higher lying Zeeman-sublevel resulting 
in a lower degree of spin polarization. In contrast the 
spin relaxation efficiency increases with increasing tem- 
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FIG. 3. (Color online) (a) Confocal laser scanning microscope image of a coffee stain ring section of sample #5. The red line 
indicates the region, which was chosen for the linear scan shown in panel (b). Here, the spectrally-integrated PL intensity is 
depicted by the solid green circles and the time-integrated DCP by the open circles. 


perature, as one can see from the shortening of the DCP 
rise time in Fig. 01(d). In the following we will first discuss 
the spin dynamics and will later focus on the dependences 
of P® q on magnetic field and NC shape. 


C. Exciton spin dynamics 

In this part the exciton spin dynamics and its depen¬ 
dence on magnetic held, temperature and NC shape is 
discussed. First, we want to stress that the observed dy¬ 
namics of the P c {t) curve of sample #1 in Fig. |4jb) orig¬ 
inates predominantly from relaxation between the dark 
exciton F = ±2 states. The plateau, which is reached at 
t 15 ns clearly indicates thermally equilibrated exci- 
tons among the Zeeman-split fine structure states. The 
characteristic rise time obtained from the fit of the ex¬ 
perimental P c (t) dependence with Eq. 0 is r s = 2.5 ns 
in this case (B = 5 T and T = 4.2 K). The thermaliza- 
tion of excitons among the bright and dark states was 
studied in detail in Ref. Notably, this study revealed 
a strong size dependence of the thermalization time that 
was found to depend linearly on the square of the NC 
radius. At temperature T = 2.6 K for NCs with a radius 
below 3 nm the thermalization was found to occur on 
timescales taf < 100 ps. Considering the range of sizes 
in our set of samples, we expect the bright to dark ther- 
malization to occur on a similar timescale. This is more 
than one order of magnitude shorter than the typical spin 
relaxation times we determine in our experiments, espe¬ 
cially at low magnetic fields. We believe that we are not 
able to resolve the dynamics on such a short timescale due 
to the limited resolution of our setup (« 250 ps). How¬ 


ever, since P c (t) at t = 0 is equal to zero, we suggest, 
that the relative population of the Zeeman-split bright 
states is equal, i.e. no thermalization takes place at such 
short timescales. As a result this equal population is 
transferred onto the dark states. Consequently, it is very 
likely that the dynamics of P c (t) at early times arises 
from the relaxation between the |+2) and |—2) dark ex¬ 
citon states. 

Figures UK a) and E^b) show the spin relaxation rate 
1/t s , obtained for sample #1 from the fit of the time- 
resolved DCP curves with Eq. 0. as function of mag¬ 
netic field and temperature. Examples of such fits are 
shown by the red lines in Figs. EKc) andQJd). As can 
be seen in Fig. [5] (a), the spin relaxation rate is strongly 
dependent on the magnetic field. At T = 4.2 K the 
rate increases from l/r s ~ 0.1 ns -1 to 3.3 ns -1 (this 
corresponds to a shortening of the relaxation time from 
t s ~ 10 ns to 0.3 ns) with increasing magnetic field from 
B = 0 T to 15 T. As one can see in the inset of Fig.EJa), 
at low magnetic fields the spin relaxation rate follows a 
cubic magnetic field dependence and it starts to deviate 
from this trend at B « 8 T. At high fields the rate seems 
to depend linearly on magnetic field. The temperature 
dependence of the spin relaxation rate is a linear func¬ 
tion, see Fig.EKb). As the temperature is increased from 
T = 2.2 K up to 12 K the spin relaxation rate increases by 
about one order of magnitude. The linear dependence on 
temperature points towards a phonon mediated process. 
For confined phonon energies much smaller than the ther¬ 
mal energy, the Bose-Einstein distribution for the phonon 
occupation can be approximated by linear temperature 
dependence: N Uq = \/(^ e hul i/ k BT _ i) ~ ksT/hujq , where 
u} q is the frequency of a phonon with wavevector q. 

In order to get deeper insight into the origin of the 
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FIG. 4. (Color online) (a) Scheme of the lowest exciton fine structure states for zero and non-zero magnetic field, (b) Typical 
polarization-resolved PL decay and resulting time-resolved DCP for sample #1. The green line is a fit according to Eq. o 
with t s = 2.5 ns. (c) Magnetic field dependence and (d) temperature dependence of the time-resolved DCP. The red curves 
are fits according to Eq. ©. 


spin relaxation we investigate its dependence on the NC 
shape. In Fig. |6ja) the magnetic field dependences of 
the spin relaxation rate are shown for all studied DiR 
samples as well as the spherical samples #R1 and #R3. 
All of them show qualitatively the same trend, i.e. at 
low fields the spin relaxation rate follows a cubic depen¬ 
dence, as shown for selected samples in Fig. Ob), while 
at higher fields a linear magnetic field dependence is ob¬ 
served. However, the spin relaxation rate seems to be 
strongly dependent on the core size and independent of 
the shell shape. This is most obvious from the green 
and blue data points in Fig. [^a), which belong to the 
spherical NCs #R1 and DiRs with a core size of 3.2 ran. 
Although their shell shapes are totally different, the mea¬ 
sured spin relaxation rates are quite similar. This ab¬ 
sence of shell shape dependence for the spin relaxation 
rate allows us to draw interesting conclusions to be dis¬ 
cussed in the following. 

The sensitivity of the spin relaxation rate on tempera¬ 
ture suggests an acoustic phonon mediated process. Due 


to the increasing number of available phonons with in¬ 
creasing temperature the spin relaxation becomes more 
efficient. However, in contrast to bulk material, the 
phonon density of states in NCs is modified by confine¬ 
ment, as it has been confirmed in many experiments ^ — 
This leads to the ” acoustic phonon bottleneck” ^ mean¬ 
ing that low-energy, long-wavelength acoustic phonons 
are suppressed in nanometer-sized structures. Due to the 
similar elastic properties of CdSe and CdS it is expected 
that the phonon confinement is different in DiRs com¬ 
pared to spherical NCs. Namely the lower confinement 
along the rod axis presumably gives rise to longitudinal 
modes in addition to the torsional and spheroidal eigen- 
modes observed in spherical structures^ This raises the 
question whether the phonon modes, which couple to 
the excitons, are considerably different in DiR structures 
compared to the spherical case. In this respect we in¬ 
fer from the data presented in Fig. jS^a) that if acous¬ 
tic phonons play a crucial role regarding the efficiency 
of the spin relaxation mechanism, the shell shape does 
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FIG. 5. (Color online) (a) Magnetic field dependence and (b) 
temperature dependence of the exciton spin relaxation rate 
1/ts for sample #1. Dashed lines are guides to the eye. 


not lead to a significantly changed number of phonon 
modes, which interact with the exciton. This conclu¬ 
sion is in agreement with the results of a recent experi¬ 
mental study employing fluorescent line narrowing (FLN) 
on DiRs in magnetic fields Therein the authors found 
that the acoustic phonon modes of DiR nanostructures 
involved in exciton recombination in presence of a mag¬ 
netic field are similar to those of confined modes in an 
isolated spherical NC. Following the preceding line of 
thought, the magnetic field dependence of the spin re¬ 
laxation rate could partly be explained in terms of the 
phonon bottleneck. For all studied samples the magnetic 
field dependence can be divided into two regimes. In 
the low field regime the relaxation rate follows a cubic 
dependence, whereas at higher fields the dependence be¬ 
comes linear. On the other hand one can claim that at 
higher fields the slope of the spin relaxation rate depen¬ 
dence is increased compared to lower fields. The onset 



Magnetic field (T) 
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FIG. 6. (Color online) (a) Magnetic field dependence of the 
spin relaxation rate l/r a for the DiR samples #l-#7 and 
spherical core/shell samples #R1 and #R3 at T = 4.2 K. (b) 
Magnetic field dependence of the spin relaxation rate 1 /t s 
for selected samples at low magnetic fields. Dashed lines are 
guides to the eye. 


of this increase occurs at B = 8 — 10 T. Assuming an 
average angle between the magnetic field and the crystal 
c-axis of 45° and taking g ex = 2.7 for the dark exciton 
g factoi'fli the corresponding Zeeman splitting amounts 
to A E 1 meV. This energy matches roughly the vi¬ 
brational energy of the l = 2 spheroidal mode of a CdSe 
sphere with a diameter of 5 nm as estimated according to 
Lamb’s theory^ We propose that the observed change 
in the magnetic field dependence results from the onset 
of these spheroidal vibrations at magnetic fields around 
8 — 10 T. Although the core diameter of the samples 
studied here is considerably smaller than 5 nm, we think 
that due to the similar elastic properties of CdSe and 
CdS, the overall size of the shell might be relevant in this 
regard^ In the Lamb theory the energy of the spheroidal 
modes scales as E oc 1 /D, where D is the sphere diam- 
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FIG. 7. (Color online) Magnetic field dependence of the spin 
relaxation rate 1 /r a for spherical NCs #R2-^R5 with differ¬ 
ent shell thicknesses and similar core size of 2.7 nm. 


eter. We therefore expect, that for larger NCs the onset 
of the linear magnetic field dependence should occur at 
lower fields. However, we are not able to verify such a 
trend from our data. Nevertheless we would like to stress, 
that our interpretation is in accordance with the find¬ 
ings of a recent magneto-optical experiment on trions in 
CdSe/CdS core/shell NCs, where the spin relaxation was 
found to be completely inhibited for small Zeeman split¬ 
ting, as a result of the acoustic phonon confinement.-^ 
The red data points in Fig. [U)a), which belong to DiR 
samples with core diameter 2.5 nm, indicate that the shell 
thickness might have a small effect on the spin relaxation 
rate. With increasing shell thickness the spin relaxation 
rate decreases slightly at a given magnetic field. In order 
to investigate this effect more systematically, we studied 
the set of spherical NCs #R2-#R5, which have similar 
core sizes but different shell thicknesses. Sample #R2 
consists of bare CdSe cores. The corresponding mag¬ 
netic field dependences of the spin relaxation rates are 
shown in Fig. [7] The most striking observation is that 
the rate of exciton spin relaxation in the bare core sam¬ 
ple #R2 is approximately one order of magnitude largeer 
(the spin relaxation time is as short as r s ft: 300 ps al¬ 
ready at low magnetic fields) compared to core/shell sam¬ 
ples (r s ft 3 ns at low fields). The addition of a CdS 
shell leads to a strong drop of the spin relaxation rate, 
but further increase of the shell thickness from 0.8 nm 
up to 5.3 nm reduces the rate only slightly. This weak 
dependence on shell thickness was also observed in a pre¬ 
vious magnetic field experiment on spherical core/shell 
ensembles^ Here, only slight changes in the spin relax¬ 
ation rate for shell thicknesses below 5 nm were observed, 
whereas for even thicker shells the decrease got more pro¬ 
nounced. This was attributed to the onset of NC charging 
at shell thicknesses around 5 nm, leading to very differ¬ 
ent spin dynamics of negatively charged excitons (trions) 



1/V (nm’ 3 ) 

FIG. 8. (Color online) Dependence of the spin relaxation rate 
1/ts on core volume V at different magnetic fields for the 
samples shown in Fig. Ha). Dashed lines are guides to the 
eye. 


compared to the neutral excitons. 

In Fig. HI we show the spin relaxation rate as a func¬ 
tion of core volume for the same samples as in Fig. (HKa). 
Interestingly, the relaxation rate seems to be inversely 
proportional to the CdSe core volume. This observation, 
along with the surprisingly fast spin relaxation for bare 
cores and the weak dependence of t s on the shell thick¬ 
ness, substantiates the importance of the NC surface and 
the CdSe/CdS interface in core/shell structures with re¬ 
spect to the spin relaxation mechanisms. We propose, 
that the spin relaxation rate depends on the strength of 
the interaction of the exciton with the surface or the in¬ 
terface. As the volume of the core increases, the surface 
to volume ratio of the core decreases and the interac¬ 
tion of excitons with the bare core surface or core/shell 
interface is continuously reduced. 

It is well established that, due to the high valence 
band offset between CdS and CdSe the hole is confined to 
the core, whereas the electron wavefunction, due to the 
smaller conduction band offset, can leak into the shell 
to a certain degreed On the one hand this could ex¬ 
plain the weak dependence of the spin relaxation rate on 
the shell thickness, since the overall surface in core/shell 
NCs is experienced only by the electron. On the other 
hand, this aspect and the strong core size dependence of 
t s indicate that the hole might be the important carrier 
regarding the exciton spin relaxation process. 

It is worth noting, that changing the core volume or 
the shell thickness of the NC always results in a modi¬ 
fication of the electron-hole exchange interaction due to 
the modified carrier wavefunction overlap. Thereby the 
exciton fine structure, and in particular the bright-dark 
energy splitting, are modified. This influences the mixing 
of bright and dark exciton states due to the perpendicular 
component of the external magnetic field. Nevertheless, 
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considering the wide range of dark-bright energy split¬ 
tings, which for the DiR samples investigated here vary 
approximately between 2 and 5 meV,— and considering 
the small change of the spin relaxation rate with shell 
thickness, we are convinced that the impact of changes 
of the exciton fine structure due to a varying electron- 
hole exchange interaction on the spin dynamics is small. 

The apparent relevance of the core/shell interface for 
the spin relaxation is interesting, because its influence on 
the optical properties of NCs is yet poorly understood but 
undoubtedly important. For example, it was shown in 
Ref. 1H, that a ” smooth transition” from the core to the 
shell material favors non-blinking emission of the NCs. 

In this context, we want to highlight the theoretical 
studies of Knipp and Reinecke as well as related work 
dealing with the interaction between charge carriers and 
acoustic phonons in semiconductor nanostructures,— 
which to the best of our knowledge have not been con¬ 
sidered explicitly in relation with colloidal NCs. In these 
studies the authors investigate a new type of electron- 
phonon coupling, which they name the ’’ripple mecha¬ 
nism”. This coupling is related to the spatial depen¬ 
dence of the electron effective mass and its potential en¬ 
ergy. The coupling occurs due to the motion of the in¬ 
terface between the nanostructure and the surrounding 
medium. Notably, the authors predict this effect to be 
dominant over the common deformation potential cou¬ 
pling for small NCs, as studied here. Moreover they ex¬ 
pect the phonon scattering rate to be strongly dependent 
on the nanostructure size with higher scattering rates for 
smaller NCs. As it was discussed for epitaxially grown 
III-V quantum dots, the ripple mechanism also strongly 
affects the electron^ and holeH spin relaxation in mag¬ 
netic field. 

The coupling between excitons and phonons in col¬ 
loidal NCs is a topic, which is still poorly understood, 
although it obviously plays an important role for energy 
and spin relaxation as well as the recombination of ex¬ 
citons. The aforementioned studies stress the relevance 
of the interface between the nanostructure and its sur¬ 
rounding medium regarding exciton-phonon coupling as 
well as spin relaxation and our experimental results point 
in the same direction. However, we think, that exciton 
spin relaxation in NCs might be result of a complicated 
interplay of many effects, since for example the presence 
of magnetic moments of defects or dangling bonds, which 
are often considered in the context of surface or interface 
related phenomena, have not been explicitly considered 
here. 


D. Equilibrium circular polarization degree 

In this section we discuss the dependence of P® q de¬ 
rived from the fit of the tinre-resolved DCP curves with 
Eq. fl3J on magnetic field, temperature and NC shape. 
In Fig. [9] the magnetic field dependence of P® q is shown 
for the spherical core/shell sample #R1 at T = 2.2 K 


and 4.2 K. At low magnetic fields P® q shows a steep in¬ 
crease, whereas at high magnetic fields it saturates. The 
blue solid lines are fits of the experimental data with 
Eq. Q, assuming a randomly orientated NC ensemble, 
i.e. foi{x) = 1. In this case, the theoretical limit for P/ q 
is expected to be —0.75, as indicated by the black dashed 
line. The fit yields an exciton g factor of g ex = 2.5, which 
is close to the recently reported values for the dark ex¬ 
citon g factor of g ex = 2.8 in a FLN experiment— and 
g ex = 2.7 in a single-NC PL experiment—, although sim¬ 
ilar values of g ex = 2.4 — 2.6 were also reported for the 
bright exciton— At a given magnetic field P/ q decreases 
with increasing temperature, see the inset of Fig. [9] An 
increase of the temperature from 2.2 K up to 25 K leads 
to a drop of |P® q | from approximately 0.5 down to 0.1. 
This is result of thermal population of the upper Zeeman- 
sublevel, corresponding to a reduced exciton spin polar¬ 
ization. The experimental data can be fitted well with 
Eq. 13]) using g ex = 2.5, see the red line. 



FIG. 9. (Color online) Magnetic field dependence of P® q for 
sample #R1 at T = 2.2 K and 4.2 K. The black dashed line 
indicates the theoretical limit of P® q = —0.75 for an ensemble 
of randomly oriented NCs. The blue solid lines are fits of 
Eq. (O to the experimental values, assuming f 0T (x) = 1. The 
inset shows the temperature dependence of P® q at B = 3 T. 
The red line is calculated taking an exciton g factor g ex = 2.5 
as determined from the fit of the magnetic field dependence. 

In Fig. flUI a) we show P/ q as a function of B for sev¬ 
eral DiR samples together with the spherical core/shell 
sample #R1 at T = 4.2 K. Apparently, the evolution 
of P/ q with B is quite different for samples of different 
shape and different core size and we always find a lower 
polarization degree for DiR NCs. The dependences differ 
in their slope at low magnetic fields and also appear to 
have very different saturation levels, although saturation 
is not fully reached within the magnetic field range of 
our experiment. A varying slope at low magnetic fields 
could result from differences of the exciton g factor, which 
is mainly determined by the variation of the electron g 
factor. In the studied structures the holes are strongly 
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FIG. 10. (Color online) (a) Magnetic field dependences of 
P® q for samples of different core size at T = 4.2 K. The black 
dashed line indicates the theoretical limit of P® q = —0.75 for 
an ensemble of randomly oriented NCs. (b) P® q for the DiR 
samples #l-#7 and the spherical core/shell sample #R1 at 
T = 4.2 K and B = 15 T as function of the aspect ratio. 
The inset schematically illustrates the corresponding change 
of NC shape. The dashed lines link data points for NCs with 
the same core sizes and are guides to the eye. 


localized in the core and the hole g factor can be esti¬ 
mated using the expression developed for a potential of 
spherical symmetry as gh = —1.09, see Refs Ha and [66l. It 
is independent of the core size and is therefore expected 
to be the same for all samples. However, the increase of 
the CdS shell leads to a decrease of the electron quantum 
confinement energy and thus to a decrease of the electron 
g factor because of its energy dependence^ Simultane¬ 
ously, the leakage of the electron wave function into the 
CdS shell leads to a contribution from the CdS g factor, 
which is close to 2, and thus to an increase of the electron 
g factor. These two effects may partly compensate each 
other. However, the latter can even result in an axial 
anisotropy of the electron g factor. 


Nevertheless, differences in the g factor do not affect 
the saturation level of P® q at high magnetic fields. Al¬ 
though saturation is not reached in the magnetic field 
range of our experiment, the data in Fig. 1101 a) indi¬ 
cate saturation of P® q at different levels. In Fig. fTOl' b) 
P® q (P = 15 T) is shown for all DiR samples as well as 
sample #R1 as function of the aspect ratio (the ratio of 
length to width of the NC shell given in Table U}. The 
values of P® q (R = 15 T) vary between —0.66 for the 
spherical sample #R1 and —0.4 for DiR sample #7. For 
a given core size, we find, that P® q is always lower for 
the NCs with higher aspect ratio. This observation indi¬ 
cates, that the differences in P® q could be related to the 
geometry of the NC shell. As discussed in Ref. j68j, the 
circular polarization of the PL in magnetic field can be 
reduced by an optical anisotropy, resulting from linear 
polarization of the individual emitters. Typical values 
for the PL linear polarization along the anisotropic c- 
axis of individual DiRs were reported to be close to 75% 
at room temperature To the best of our knowledge, 
there are no reports about the linear polarization for sin¬ 
gle DiR NCs at low temperatures. The high linear po¬ 
larization degree of the PL observed for single DiR NCs 
at room temperature may result partly from the exciton 
fine structure and partly from a dielectric enhancement 
effect ; 20 ’ 69 both of which are dependent on the NC ge¬ 
ometry. In the following section we provide a detailed 
discussion of these effects and analyze how they affect 
the circular and linear polarization of the PL in magnetic 
field for NCs with various shapes at low temperatures. 


IV. THEORETICAL TREATMENT OF LINEAR 
AND CIRCULAR POLARIZATION OF LOW 
TEMPERATURE PL IN EXTERNAL MAGNETIC 
FIELD 

A. Model 

We start with the assumption, that the lowest exci¬ 
ton state is the dark exciton with momentum projection 
|P| = 2 on the hexagonal c-axis (which is also the rod 
axis in DiR NCs). Figure fill shows a sketch of the ex¬ 
citon energy levels in zero and nonzero magnetic Held. 
Only the ground state F = ±2 and the bright exciton 
states F = ±1 L ' U and F = 0 U are shown (without the 
energy splitting between the l L and l u states). Fur¬ 
thermore, we assume that the dark exciton IF 1 ) = 2 can 
recombine radiatively only via coupling with the bright 
exciton states. Coupling with the \F\ = 1 L - U states leads 
to recombination which polarization is oriented isotropi¬ 
cally in a plane perpendicular to c-axis^ The respective 
recombination rate is denoted by r 0 ! at zero and r 2 i at 
nonzero magnetic field so that the probability to emit 
light in the direction k (see the schematic of the exper¬ 
imental geometry Fig. 12) is given by r° 1 (l + cos 2 0) 
and r 2 i(l +cos 2 0), respectively. Moreover the coupling 
with the \F\ = 1 states can be enhanced by a magnetic 
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FIG. 11. (Color online) Schematic of the band edge exci- 
ton energy levels in zero and nonzero magnetic fields. Only 
the dark exciton ground state with F = ±2 and the bright 
exciton states F = ±1 L,U and F = 0 U are shown (without 
splitting between the 1 L and l r states). The red and blue 
arrows show the coupling pathways, which may enable dark 
exciton radiative recombination. The sketches at the bottom 
illustrate the respective orientation of the transition dipole for 
the different coupling pathways with respect to the crystal c- 
axis. Here colors correspond to those of the arrows indicating 
the coupling pathways. 


field oriented perpendicular to the c-axis j 15 ’ 44 For this 
reason, r 2 i depends on the angle 0 between the mag¬ 
netic field direction and the c-axis (see Fig. [12]). In con¬ 
trast, the coupling of the dark exciton |F| = 2 with the 
bright F = 0 U state leads to recombination polarized 
parallel to the c-axis. Its recombination rate, which is 
denoted by r 2 o, is independent of the magnetic field and 
the probability to emit light in the direction k is given 
by r 2 osin 2 0. The total recombination rate of the dark 
exciton averaged over all light directions is proportional 
to r 20 + 2r 2 i. 

We consider hereafter only the range of magnetic field 
satisfying the condition g e xFBB sat <C A Em- In this 
range the recombination rates of F = —2 and F = +2 
excitons calculated in the framework of second order per¬ 
turbation theory are equal. 

It should be noted that the usual assumption leading 
to Eq. J3]) for describing the magnetic field dependence of 
P® 9 is that the dark exciton recombines due to coupling 
with the |F| = 1 L,U bright state s 15 ! 44 , characterized by 
r 2 i and thus has the same polarization properties as the 
|F| = 1 exciton. In this case, the dipole moment ac¬ 
quired by the dark exciton is oriented transverse to the 
c-axis, see the scheme left bottom in Fig. [TT] The new as¬ 
pect of our consideration is the assumption, that the dark 
exciton can recombine with the properties of the F = 0 
bright exciton as well. This coupling between the |F| = 2 
dark exciton and the F = 0 state described by r 2 o could 
be mediated by l = 2,m = ±2 acoustic phonon s 71 ' 72 . It 


FIG. 12. (Color online) Schematic of the experimental ge¬ 
ometry for measuring circular and linear polarization of the 
PL from DiRs or spherical NCs with anisotropic c-axis. Note 
that in the ensemble the angles 0 and <3>, which determine 
the direction of c with respect to the light vector k and the 
linear polarizer axis ey, are arbitrary. 

could be also caused by the shape anisotropy in the plane 
perpendicular to the c-axis, which mixes the valence band 
states^ but more importantly leads to an anisotropy of 
the electron wavefunction leaking into the rod and thus 
to an anisotropy of the Coulomb interaction in the exci¬ 
ton. 

Importantly, in DiR NCs the recombination polarized 
along the c-axis can be additionally enhanced by the di¬ 
electric enhancement effect. The enhancement factor R e 
in Ref. H3 is given by 

Km + (k s - K m )n ( ' x ’ y ' > \ 2 
Km + [k s - K m )n (*) ) 

where K m and k s are the dielectric constants of the sur¬ 
rounding medium and the semiconductor, respectively. 
n O,i/) = (x _ 7 j( 2 ))/2 and = a 2 [ln(26/a) — 1 ]/b 2 are 
the depolarization coefficients along the directions per¬ 
pendicular (x,y) and along (z) the hexagonal axis. The 
NCs are approximated by ellipsoids with major semi-axis 
b and minor semi-axis a. Essentially, the rod shaped shell 
of DiRs acts as a linear polarizer, which affects the ab¬ 
sorption as well as the emission of light. Note that here 
we neglect the frequency dependence of R e and assume 
that it is the same for excitation and emission. 

The dielectric enhancement effect can lead to selective 
excitation of DiR NCs oriented perpendicular to the mag¬ 
netic field (the light propagation direction). According to 
Ref. [73, the probability for a NC to be excited by linearly 
polarized light is 1 + ( R e — l)(e ■ c) 2 , where e is the light 
polarization vector. For the geometry shown in Fig. [T2] 
and light polarized along the e\\ direction one can write 
(e • c) 2 = sin 2 0cos 2 <1>. By averaging over all polariza¬ 
tions in the plane (e||,ej_) we obtain the probability for 
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excitation of a NC, oriented at an angle 0 with respect 
to the light propagation direction, by unpolarized light 
as 


jexc(x) = 1 + 0.5(P e — 1)(1 — x 2 ) , (6) 


where x = cos@. Taking into account the selective ex¬ 
citation effect, the general expression Eq. ©, which de¬ 
scribes the DCP for an ensemble of DiR NCs, transforms 
into 

p _ Jo [Irr+{ X ) ~ { x )]foi( x )f exc (x)dx 
J^[I a+ (x) - I a -{x)]f OT {x)f exc {x)dx 

In the following we consider only the case of randomly 
oriented NCs, i.e. / or (x) = 1. When the selective excita¬ 
tion is taken into account and the I^o rate is neglected, 
expression Eq. (0 for the magnetic field dependence of 
the DCP remains unchanged with the function / or (x) re¬ 
placed by /excfz)- 

When we additionally consider the activation of the 
dark exciton by a magnetic field applied perpendicular 
to the c-axis, which leads to a field and angle dependence 
of the rate T 2 i(B,x), Eq. 0 transforms into 

^ _ Jq 1 'jjB, x)2x ta,nh(AE/2kBT)f exc (x)dx 

fy ry(B,x)(l + x 2 )f exc (x)dx 

( 8 ) 

where the factor y(P, x) = P 2 i (P, x^T^ is a measure of 
the strength of the magnetic-field-induced mixing. 

In general, the angular and magnetic field dependence 
of T 2 i(B,x) can vary, depending on the ratio between 
g ex PBB sat and the bright-dark exciton splitting A Ew 
B sat is the magnetic field at which the DCP is saturated. 
For the range of magnetic fields we consider here the 
condition g ex PBB sat -C A Em is satisfied and one can 
write T 2 i(x) = r ^![1 +p(l — x 2 )], where the parameter p 
describes the dependence on B with p oc B 2 j 15 i 75 

Now we also include the coupling of the \F\ = 2 states 
to the state F = 0 U . The transitions polarized along 
the c-axis contribute equally to I a + (x) and I cr _ (x) with 
the probability r 2 oP e sin 2 0. Accounting for these tran¬ 
sitions we obtain the following expression for an ensemble 
of randomly oriented DiR NCs: 


^ 'y(B,x)2xta,nh(AE/2k B T)f exc (x)dx 

Jo b/( B > x )0- + x 2 ) + r(l - x 2 )\ f exc (x)dx ’ 

(9) 

where r = P e r 2 o/r 2 i- Obviously, on the one hand the 
coupling to the F = 0 U state lowers P® q at a given mag¬ 
netic field. On the other hand it leads to a higher degree 
of linear polarization for the individual NC, as will be 
discussed below. 

The degree of linear polarization for an individual NC, 
is pi = (I || — I±)/(I\\ + / 1 ). Here T||(Tj_) is the intensity 
of light polarized along the e||(ej_) direction. It can be 
shown that 


(i - x 2 )(f? e r 20 - r 2 i) cos( 2 $) 
Pl[X ’ ’ (l + x 2 )T 21 + (l~x 2 )R e r 20 


( 10 ) 


One sees that after averaging over the randomly ori¬ 
ented ensemble of DiR NCs, the overall linear polariza¬ 
tion Pi = f^dx pi(x, 4>)d<h is zero. However, it still 
can be nonzero for the individual DiR NC reaching a 
maximum for rods oriented at d* = 0° and 0 = 90° (i.e. 
x = 0). In the following we consider only the maximum 
linear polarization for this particular geometry which for 
zero magnetic field gives 


Pi,o' 


r — 1 
r + 1 


( 11 ) 


It is worth noting, that for the case r < 1 (i.e. a higher 
coupling rate to the ±1 states T 21 in comparison to the 
dielectrically enhanced rate P e P 2 o), pTo* has negative 
values. This corresponds to the case, where the light is 
preferentially polarized perpendicular to the c-axis. In 
magnetic fields at which the circular polarization degree 
is saturated (B = B sat ): 


„max ( td sat \ 

Pl ) 


r — (1 +p sat ) 
r + (1 + p sat ) 


( 12 ) 


Here we have defined (1 + p sat ) = 7 (P sat ,x = 0) = 
r 2 i(P sat ,x = 0 )/r 21 • It is a measure for the increase 
of r 2 i at B = B sat for NCs oriented perpendicular to the 
magnetic field. 

So far, the equations for the circular polarization of 
the DiR ensemble and the linear polarization of individ¬ 
ual DiRs in magnetic field have been derived. We turn 
now to the analysis of the correlation between the linear 
polarization of individual DiRs and the ensemble DCP in 
more detail and compare our theory to the experimental 
data on DiRs. 


B. Discussion 

In the following, we analyze how the saturation level 
of P® q (i.e. P aat ) is affected by the different mechanisms 
we introduced in the previous section. To that end, we 
simplify the discussion by analyzing first the effect of 
selective excitation and the newly introduced coupling to 
the F = Q u state, while we set p = p sat = 0, which means 
that magnetic field activation of dark exciton radiative 
recombination is neglected as first approximation. Later 
on the additional effect of p ^ 0 is discussed. 

In Fig. Q0a) we show the magnetic field dependences 
of P® q for the DiR samples #l-#3 and the spherical 
core/shell sample #R1 at T = 4.2 K. Obviously, the 
dependences are quite different and the P® q tend to sat¬ 
urate at different values. The lines are calculations ac¬ 
cording to Eq. 0. The parameters used are g ex = 2.7, 
R e = 3.2,3.3,3.4,1 and r 20 /r^ = 0.03,0.13,0.19,0.02 
for samples #1,#2,#3 and #R1, respectively. This 
means that we assume the ratio between the coupling of 
the dark exciton to the F = 0 state and the coupling to 
the F = ±1 states to be different for the different NCs. 
The enhancement factors R e were calculated according 
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FIG. 13. (Color online) (a) Magnetic field dependences of P aq for the DiR samples #1, #2, #3 and the spherical core/shcll 
sample #R1 at T = 4.2 K. The dashed and dotted lines are calculated according to Eq. ([9]). (b) Saturation values of the 
circular polarization in high magnetic fields P® at for an ensemble of arbitrarily oriented NCs vs. maximal linear polarization of 
individual NCs in zero magnetic field p™Q X . The orientation of the transition dipoles with respect to the NC c-axis is illustrated 
in the sketches below. The black curve is calculated without considering the effect of selective excitation (i.e. R e = 1) and the 
red curve is calculated taking R e = 1000, which corresponds to an extremely strong effect of selective excitation. The inset 
shows P c sat a s a function of R e for p™= — 1 (blue curve) and pTo* = 0 (green curve). The dashed lines in the inset mark the 
range of calculated enhancement factors for the investigated samples according to Eq. ©. The dashed and dotted lines in the 
main plot of panel (b) are the saturation values derived from the corresponding calculated curves in panel (a). All calculations 
were done neglecting the effect of magnetic field activation of the dark exciton recombination. 


to Eq. ©, where the dielectric constants where taken to 
be n m = 2 and k s = 6 for the surrounding medium and 
the NC, respectively^ Apparently, the experimentally 
measured P® q can be well reproduced by our model. 

Figure HHb) shows a plot of the saturated ensemble 
polarization P® at (for randomly oriented DiRs) as func¬ 
tion of the maximal linear polarization of individual DiRs 
at zero magnetic field The magnetic field activa¬ 

tion of the dark exciton is still neglected (i.e. p sat = 0). 
The horizontal lines correspond to the saturation values 
of P® q derived from the corresponding calculated curves 
in the left panel (the black arrows indicate saturation). 
The red and black curves are the theoretically predicted 
dependences of P| at on p™ q x with and without consider¬ 
ing the effect of selective excitation. Obviously, the se¬ 
lective excitation effect leads to an additional decrease of 
the saturated polarization value. The dependence of P aat 
on R e due to selective excitation is shown in the inset for 
p] r \f x = — 1 and pfjf* = 0 (see the corresponding blue and 
green arrows). The dashed lines indicate the range of en¬ 
hancement factors calculated for the DiRs investigated 
here. Notably, even in the case of an extremely strong 


dielectric enhancement effect, the saturation value is not 
lowered by more than approximately 12%. In contrast, 
as can be seen from the red and black curves, the increase 
of the r value (i.e. the ratio of the coupling rates, see the 
upper scale in Fig. 113T b)) and the resulting maximal lin¬ 
ear polarization of individual DiRs at zero magnetic field 
leads to a decrease of the saturation value of the circu¬ 
lar polarization in the ensemble that can vary between 
—0.75 and 0. From the intersection of the dashed and 
dotted lines (the theoretically derived saturation values) 
with the black curve, we can predict within this approxi¬ 
mation the maximal linear polarization of the individual 
rod at zero field to be negative for all DiRs considered 
here (Pi\° /X < 0 ). This means that we expect the light 
emitted from these DiRs to be linearly polarized predom¬ 
inantly transverse to the c-axis at low temperatures, see 
the sketches below Fig. fl3lb). 

Next, we discuss in more detail how the magnetic field 
activation of the dark exciton influences the polarization 
properties of the DiR emission. First, we analyze the 
effect on the linear polarization of individual rods. In 
Fig. M we plot the maximal linear polarization of indi¬ 
vidual rods at B = B sat versus the maximal linear polar- 
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FIG. 14. (Color online) Maximal values of individual NC 
linear polarization in high magnetic field p™ ax vs. maximal 
values of the linear polarization in zero magnetic field p“o' x 
for p sat =0,1,2,3. 


ization at zero magnetic field for different values of p sat . 
We remind, that the parameter p was assumed to be de¬ 
pendent on the magnetic field. Since we consider here 
the linear polarization either at B = 0 or B = P sat , we 
use p(B sat ) = p sat = 7(B sat ,0) — 1 as model parameter. 
It is defined by the coupling of the F = ±2 states to the 
F = ±1 states at zero and high magnetic fields and for a 
NC with the hexagonal axis oriented perpendicular to the 
magnetic field. The magnetic field enhances the coupling 
to the F = ±1 states and therefore the amount of tran¬ 
sitions polarized transverse to the c-axis increases with 
increasing magnetic field. This effect can be seen from 
the curves shown in Fig. [Td] If p sat = 0 the maximal lin¬ 
ear polarization is independent of the magnetic field, i.e. 
p™ ax = p™ ax . However, with increasing p sat the linear 
polarization in high magnetic fields shifts continuously 
to negative values, which means that the transitions po¬ 
larized transverse to the c-axis are enhanced. This can 
lead to a situation in which at zero magnetic field the PL 
is polarized parallel to the c-axis, whereas at high fields 
it is polarized transverse to it. 

The parameter p sat is unknown for the samples inves¬ 
tigated here but could be determined by measuring the 
linear polarization of an individual DiR at zero and high 
magnetic fields. Such an experiment goes beyond the 
scope of this study. However, we can estimate the or¬ 
der of magnitude of p sat by considering the shortening of 
the PL decay typically observed in magnetic field exper¬ 
iments. In Ref. 17, which reports on magneto-optical PL 
measurements on single spherical CdSe/ZnS core/shell 
NCs, the authors observed a shortening of the PL decay 
from approximately 180 ns to 50 ns for a NC oriented at 


75° with respect to the magnetic field. This corresponds 
to a ratio of 3.6 for the recombination rates and therefore 
p is equal to 2.6. We expect a slightly larger shortening 
for larger orientation angles, thus the upper limit of p sat 
should be around 3. It is worth noting, that this ap¬ 
proximation works only, if one assumes, that there is no 
contribution from the F = 0 U states to the PL, which can 
reasonably be assumed for NCs, where the DCP satura¬ 
tion value in a randomly oriented ensemble approaches 
—0.75 at high magnetic fields. 

Now, we analyze how the increased transverse linear 
polarization due to magnetic field activation of dark exci- 
ton radiative recombination affects the circular polariza¬ 
tion of the ensemble PL. In Fig.[T5]we plot the saturation 
value of the ensemble circular polarization at B = B sat 
as function of the individual NC linear polarization at 
B = 0 for p sat = 0 (black curve) and p sat = 3 (green 
curve). The selective excitation effect is not considered 
here (i.e. R e = 1). The horizontal lines are similar to 
those in Fig. fT3lb) and correspond to the predicted sat¬ 
uration levels for the samples shown in Fig. [T3l a). They 
are given here for better comparability. It is worth not¬ 
ing, that with accounting for the angular dependence of 
the dark exciton recombination rate, the predicted level 
of the saturateion DCP values from fitting the magnetic 
field dependences P® q (P) will be lower for each sample. 
One can see that in the region of small values of r (see the 
upper scale of Fig. USD, the effect of the dark exciton mag¬ 
netic field activation leads to a decrease of P® at averaged 
over the ensemble. This is because the magnetic field 
activation depends on the angle and the resulting tran¬ 
sitions become stronger for the NCs oriented perpendic¬ 
ular to the c-axis where the Zeeman splitting is smaller. 
With the increase of r, however, the increase of the T 2 i 
transition becomes more important than the angular de¬ 
pendence and the resulting P® at is increased. As a result, 
from the intersection of the dashed and dotted lines with 
the green curve for p sat = 3, we can predict the maximal 
linear polarization of an individual rod at zero magnetic 
field. It is positive for the DiRs #2 (p™g x = 0.10) and 
#3 (p™q x = 0.32), so that we expect the light emitted 
from these DiRs to be linearly polarized predominantly 
along to the c-axis at low temperatures. However, as can 
be seen from Fig. 1141 the maximal linear polarization of 
an individual rod at magnetic field B = B sat is negative: 
p™ = —0.53 and —0.38 for the samples #2 and #3, 
respectively. 

As shown above, our model, which introduces the cou¬ 
pling of the dark exciton to the bright exciton state 
F = 0 r/ in addition to the well-known coupling to the ±1 
states, allows one to explain why for DiR NCs a signifi¬ 
cantly lower DCP is observed in magnetic fields compared 
to spherical NCs. The observed dependence on the aspect 
ratio highlights the role of the dielectric enhancement ef¬ 
fect, which favors emission polarized parallel to the c-axis 
in DiRs. The reason for the additional increase of the ra¬ 
tio r 20 with increase of the DiR aspect ratio can be 
explained by the increase of the shape anisotropy in the 
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FIG. 15. (Color online) Saturation values of the circular po¬ 
larization degree in high magnetic fields P/ at for an ensem¬ 
ble of arbitrarily oriented NCs vs maximal linear polarization 
values of individual NCs in high magnetic fields p™Q X , tak¬ 
ing into account the effect of magnetic field activation of the 
dark exciton recombination (p sa = 0 for the black curve and 
p sat = 3 for the green curve). Calculations are done without 
considering the effect of selective excitation ( R e = 1). 


plane perpendicular to the rod axis or by the increase 
of the Coulomb interaction anisotropy. Admittedly, the 
suggested model has a few free parameters and the the¬ 
oretical prediction of the linear and circular polarization 
properties for DiR have to be considered with care. Ad¬ 
ditional dedicated measurements at low temperatures are 
necessary to elucidate the relevance and the magnitude 
of each parameter and their dependences on NC size and 
shape. Especially the determination of the linear polar¬ 
ization properties of individual DiRs at low temperatures 
and in magnetic fields would be highly enlightening in 
this regard. Nevertheless, we believe, that these find¬ 
ings may lead to a deeper understanding of the processes 
related to the exciton fine structure and the resulting 
polarization properties of the emitted light. Finally, we 
would like to point out, that also for spherical core/shell 
NCs saturation DCP values considerably lower than the 
theoretically expected —0.75 have been observed in high 
magnetic field experiments and these low DCPs could 
not be explained at that time i 76 ' 77 However, they could 
be explained in the framework of our theory, if one con¬ 
siders, that also in spherical core/shell NCs the coupling 
to the F = 0 r/ state is relevant. This would intrinsically 
lower the saturation value of the DCP in the ensemble as 
discussed in connection with Fig. 1131 Furthermore, real 
NCs are never perfectly spherical, which makes it nec¬ 
essary to consider also the dielectric enhancement effect 
for these structures. 


In summary, we performed a comprehensive study 
of the PL polarization and exciton spin dynamics in 
CdSe/CdS core/shell NCs in external magnetic field. We 
demonstrated the possibility of probing the orientation 
of NCs with anisotropic shell shape with respect to the 
magnetic field axis by measuring the PL DCP induced 
by the field, and highlighted possible applications of this 
technique for testing nanorod based arrays. 

Regarding exciton spin dynamics we found that the 
spin relaxation of dark excitons is strongly related to the 
surface of CdSe in bare NCs or the interface between 
CdSe and CdS in core/shell systems. This was evidenced 
by a large decrease of the spin relaxation rate after adding 
a CdS shell. The importance of the interface in core/shell 
systems was highlighted by a strong core size dependence 
of the spin relaxation rate, which decreases for increasing 
core volume. In contrast, the overall NC surface seems 
to play only a minor role in this regard, as concluded 
from the weak dependence of the spin relaxation rate on 
shell thickness. The temperature and magnetic field de¬ 
pendences of the spin relaxation rate emphasize the im¬ 
portance of confined acoustic phonons in the relaxation 
and recombination of dark excitons. In this respect, DiR 
NCs show a similar behaviour as spherical core/shell sys¬ 
tems. We believe, that these findings may stimulate fur¬ 
ther discussion about the role of the core/shell interface 
for the optical properties, in particular the dark exciton 
recombination and spin relaxation in NCs. With respect 
to applications our results on spin dynamics may help 
to figure out the relevant parameters for spin-preserving 
environments in spintronics devices, e.g. for spin storage. 
In respect thereof we can state, that adding a shell for 
surface passivation and the use of large core NCs as well 
as thick shells favor long spin lifetimes. 

Finally, we extensively studied the circular PL polar¬ 
ization of DiRs in magnetic field. Despite a remarkably 
low equilibrium DCP for randomly oriented DiRs com¬ 
pared to spherical core/shell NCs our results point to¬ 
wards the dependence of the saturated DCP on the NC 
aspect ratio. We elaborated a theory for analyzing and 
predicting the circular and linear polarization properties 
of DiR emission. It includes the dielectric enhancement 
effect as well as dark exciton recombination via coupling 
to bright exciton states with radiative transitions polar¬ 
ized parallel or perpendicular to the hexagonal c-axis. 
This model can explain the strongly varying (almost sat¬ 
urated) DCPs typically observed even for spherical NCs 
in high magnetic fields and may lead to a better under¬ 
standing of the exciton fine structure and its relevance 
for the PL polarization properties of NCs especially at 
low temperatures. 
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